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Abstract 



We generalise our previous observation on the invariance of the Jarlskog 
determinant to matter effects in neutrino oscillations. Within the context of 
standard neutrino oscillation theory with matter effects, we present the complete 
set of (five) matter invariant observables for neutrino propagation in matter. We 
give some examples of their application. 
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1 Introduction 

Neutrino oscillations violate lepton flavour conservation. The amplitude A a p for a 
neutrino of flavour a to be detected as a neutrino of flavour (3 is given as a function 
of propagation distance L by the (matrix) equation: 

A = exp(-iHL) (1) 

where H is the effective Hamiltonian in a flavour basis (ie. in a weak basis in which the 
charged lepton mass matrices are diagonal - this is defined up to an arbitrary phase 
redefinition of the charged lepton fields). We may take: 

H=Udiag(p,A 21 ,A 3l )rf (2) 

where U is the conventional MNS matrix and the A^ are defined by Ay = (mf — m3)/2E, 
with rrii the neutrino masses and E the neutrino energy. In this paper, all calculations 
are assumed to be in such a flavour basis. 

We will find the following theorem useful. Any function, f(H), of a Hermitian 
operator H, can be expanded using Lagrange's formula [|l|]: 

f(H) = J2f(X i )X i (3) 

i 

where the A, are the eigenvalues of H and the X 1 are hermitian (projection) operators 
given by: 

X i = YltfijH- Aj) ^ 
Uj^i(\ - A,) 

Since H satisfies its own characteristic equation, the complete product Yli(H — Aj) 
(with no factors dropped) clearly vanishes. It follows that the X % are both left and 
right eigenvectors of H with eigenvalues Aj, ie. (H—Xi)X l = X l (H — Xi) = 0. With the 
normalisation defined as in Eq. |], we have {X 1 ) 2 = X 1 , so the X 1 are indeed projection 
operators, representing the (pure) mass-eigenstate density-matrices. 

The ancient formula Eq. ^ has been used previously for the specific case of three- 
family neutrino oscillations in eg. Ref. 0, and has recently re-appeared in the same 
context in Ref. 0. Applied to Eq. [I] it gives: 

A = ^X i exp(-iXiL). (5) 

i 

Any X % is readily computed from the corresponding column of the MNS matrix 0: 

Xifi = UaiUZi (6) 

(where no summation over % is implied) . We note here that the elements of the matrices 
X 1 are dependent on the phase-convention chosen for the charged lepton fields and are 
therefore not observables (likewise, of course, neither are the MNS matrix elements 
themselves) . 
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Using Eq. |3] we may compute arbitrary powers of H as follows: 

H'=J2KX\ (7) 

i 

Given the eigenvalue spectrum, Eq. ^ determines the matrices X % in terms of H, whilst 
setting r = 1 in Eq. [7| yields the inverse relation, determining H in terms of the X % . 
The case r = of Eq. |7| expresses the unitarity of U. For an n x n Hamiltonian, 
there are only n independent equations of the form Eq. [7], because the characteristic 
equation can be used to write higher powers r > n in terms of the n lowest powers 
H r , < r < n. 

The squared amplitude |A Q/ g| 2 for a ^ [3 (a = (3) gives the appearance (survival) 
probabilities as a function of propagation distance (the dependence on neutrino energy 
is implicit in the factors Ajj): 

P{a^ (3) = \A aP \ 2 = 5 Q/3 -4^^sin 2 (A y -L/2) 

+ 8J a/ ?sin(A 12 L/2)sin(A 2 3L/2)sin(A3ii:/2) (8) 

where 

K% = Re(XVQ), (9) 

and 

J aP = Im(Z^X^) = ± J (for a ^ 0) or = (for a = /?). (10) 

J is Jarlskog's CP- and T- violating invariant 0. Inspection of Eqs. |8|-p~0| indicates that 
the phase-convention- independent products of the form X^X 1 *^ {i ^ j) are observ- 
ables of the system 0, the real parts K^o parameterising CP-conserving oscillation 
amplitudes and the imaginary parts J a/ g parameterising CP violation. We note that 
unitarity imposes many constraints among the observables, K^* (and J), and that 
once three of these are specified (such that at least two different flavour indices are 
used) along with J, the remaining ones may be determined via unitarity. 

When neutrinos propagate in matter, the Hamiltonian in the flavour basis is mod- 
ified: H -> H', with 

H' = H + diag(a, 0, 0) (11) 

where a = \f2GpN e . This modifies both the eigenvalues and eigenstates of the neutrino 
Hamiltonian (in this paper, "primed" symbols will be used to denote the values of 
quantities when neutrinos propagate in matter and unprimed symbols for their vacuum 
counterparts) . 

We recently pointed out that the determinant of the commutator of the lepton 
mass matrices is invariant with respect to matter effects, 

Det[£> 2 , H'} = Det[L> 2 , H] (12) 

where is the charged-lepton mass matrix, which is diagonal by construction in this 
basis. This allowed the derivation of a matter-invariant quantity, expressible both 
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as a function of elements of the effective Hamiltonian, and as a combination of the 
eigenvalue differences, A^, and the CP violation parameter, J |7]]: 

lm{H' eil H'^H' Te ) = lm(H e ^ T H Te ) 

= A' 12 A 23 A' 31 J' = A 12 A 23 A 3 iJ. (13) 

This quantity has found application in several instances |§. It is not however, the 
only matter-invariant combination of neutrino oscillation observables. Here we count 
the number of independent such combinations and derive a complete set. 

2 A Complete Set of Matter Invariant Observables 

Clearly, in a flavour basis, the modification to the effective Hamiltonian induced by 
matter, Eq. |Tl|, affects only the "ee" element: 

H' ee = H ee + a (14) 

whereas the remaining eight elements of H are matter-invariant in this basis, ie. 

H' a(3 = H aP , a/?^"ee". (15) 

The invariants of Eq. |1| are not immediately useful, however, as the Hamiltonian 
(vacuum or matter) is not itself observable. It transforms under (separate) phase 
redefinitions of the flavour eigenstates: 

H -> (16) 

where 

$ = exp{z diag(0 e , <f>^ <p T )} (17) 

so that 

H aP - e^-^H af3 (18) 

without any observable consequence. The X 1 transform similarly. If our invariants are 
to be readily applicable, we must seek to relate observable quantities directly. 

We can count the number of independent observables represented by H as follows. 
An arbitrary n x n Hermitian matrix comprises n 2 independent real parameters (n 
on-diagonal, n(n — 1), off- diagonal). Eqs. [I~6| - |T8| imply that for the effective Hamil- 
tonian, n — 1 of these are unobservable relative phases, so that we would expect, at 
most, n 2 — n + 1 observable quantities. In addition however, subtraction of any mul- 
tiple of the identity from H (a "trace transformation") leaves all neutrino oscillation 
observables invariant, so that only eigenvalue differences are relevant, leaving finally 
only n(n — 1) independent observables in neutrino oscillations. For the 3x3 MNS 
oscillation phenomenology, we have six observables, which are conventionally consid- 
ered to be three real angles, one phase and two mass-squared differences. They could 
equally well be taken to be three independent K% (see Eq. ||), J and two eigenvalue 
differences, A21 and A31, or alternatively six phase- and trace-invariant functions of 
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elements of H (see below). When neutrinos propagate in matter, one parameter is 
modified, Eq. [14|, so that we may expect, in general, n{n — 1) — 1 matter invariances 
(eg. five, in the three-family case). 



We seek relations analogous to Eqs. PH15| but involving only functions of the 
H a p which are invariant under both phase transformations and trace transformations. 
Taking Eq. [13] as a model, we seek furthermore to write them in terms of the observables 
of Eq. || Six independent such relations should exist. Relations which do not involve 
the "ee" flavour label will, in addition, be matter- invariant. The off-diagonal elements, 
H a p (ot (3), are already invariant under trace transformations, and we can form 
phase-transformation invariants simply by taking their moduli-squared: 

\K/ = \H aP \ 2 (19) 

= \J2Kx':p\ 2 = izvc/ (20) 

i i 

= -E4^ = ~J2 A %K% a (21) 

j<k j<k 

where we have used the relation, Eq. ^ with r = 1 to obtain Eq. ^0], and with r = 
(unitarity) to obtain Eq. |21]. The quantities in Eq. ^1] are all observable (see eg. Eq. §) 
and the three independent combinations are each matter-invariant (only three combi- 
nations because \H a p\ = \Hp a \). [] 

In a complementary fashion, the diagonal elements of H are already phase-invariant 
and we can form trace-transformation invariants by taking differences, whereby: 

— H TT = — H TT 

which is matter- invariant, and 

H' ee -Tr(H')/3 = H ee -Tr(H)/3 + 2a/3 
= (X' e l - i)A' 21 + (X 3 e - i)A' 31 = (Xi - i)A 21 + {Xl - ±)A 3 i + 2a/3 (23) 

which may be said to be matter-covariant. 

The three invariances of Eq. |20] correspond to the "sum-rules" of Ref . |§ , although 
they are derived in a different way here, without appealing to the commutator of the 
mass-matrices Q. Their expression in terms of our observables, Eq. |2l|, and the relations 
of Eqs. and ^3] are new. Our derivation demonstrates that all the above invariances 
follow straightforwardly from the fact that the Wolfenstein term affects only the "ee" 
element of the effective Hamiltonian, Eqs. |I3L |T5| . 

1 The diagonal analogue of Eqs. |l9r|2l| : - Ylj<k ^jt^oa = (H aa ) 2 — (H 2 ) aa is not independent, 

A 2 rs-jk , y> a 2 K ]k 
■>3<k ^jk^-aP ' 2^,j<k ^jk^-a-fi 



the LHS equating through unitarity to: V. <fe A 2 K J a k + V <fc A 2 K^ a ^ /3 ^ 7, and the RHS 



being simply ~\H Q p\ 2 ~ \H aj \ 2 

2 The quantities Zij and Zij defined in Ref. |9) are in fact simply the off-diagonal elements of the 
effective Hamiltonian in vacuum and in matter respectively. 
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The five combinations Eqs. |2T|-|2"3| all relate CP- conserving observables (the diagonal 



elements X\ 



| 2 are observable). The combination Eq. [13| parameterises CP 
violation and is therefore independent. It is manifestly matter-, phase- and trace- 
transformation invariant, as required. We have therefore completed the derivation 
of six independent identities, Eqs. and ^T| - |23| , relating particular combinations of 
matter-modified neutrino oscillation observables to their vacuum counterparts. Five 
of the combinations (Eqs. [13], and |22[) are matter- invariant, while one combination 
(Eq. |23p depends in a simple way on the matter-density. As regards neutrino oscillation 
phenomenology, these six combinations constitute a complete set, whereby any further 
matter- invariants (which can certainly be found, eg. R,e(H efl H^ T H Te )) must always be 
expressible in terms of those given here. 



3 Some Applications 

As a first application, we show briefly how the formalism presented in the last Section 
applies to the simple case of two neutrino families. In this case, there is only one 
independent off-diagonal element of the effective Hamiltonian, \H efl \, which leads to 
the analogue of the matter-invariant in Eq. |21| : 

I ttI 1 2 I TT 1 2 

\ rL e^\ ~ \ S2 eiJ,\ 

= -A' 2 < 2 = -Aj 2 KH 
=^ A /2 2 sin 2 20' = A 2 2 sin 2 2# (24) 

while the two-family analogue of Eq. |23| is 

«e-^)A' 21 = a/2 + (A 2 e - I)A al 

^(sinV-i)A' 21 = a/2 + (sin 2 #-i)A 21 

=^-cos2fl'A 21 = a-cos2#A 21 (25) 

where we have used the definitions of the X 1 and K%g in terms of the 2x2 mixing matrix 
elements to recover the standard results We have found one matter- invariant 

combination of observables [T1J, and one which is matter-covariant, the complete set 
for the two-family case (compared with five and one respectively for the three family 
case). 

In the case of three families, we consider the appearance probabilities Eq. ||] with 
a 7^ (3, in matter of arbitrary (uniform) density. Expanding each term in A'^L/2 for 
small, we find: 



P'(a-/3) = -E<S( A ^ 2 - A ^ 4 /24 + ...) 



Kj 



+ J' a/3 A' n A 23 A' 31 L* + ... 

= |#;/L 2 + lm(H' e ^ T H' Te )L s + 0{A%L') + ... 

= \H a p\ 2 L 2 + lm(H e ^H^ T H Te )L ;i + 0(A^L 4 ) + ... (26) 
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We have used Eqs. ^T] and [13] to identify the coefficients respectively of the CP- 
conserving term of order L 2 and the CP- violating term of order L 3 with manifestly 
matter-invariant combinations of elements of the effective Hamiltonian. Thus, they 
are determined directly in terms of elements of the vacuum Hamiltonian. The term of 
order L 4 is not matter- invariant, and is the lowest-order term through which matter 
effects enter. The well-known || result that matter has no observable effect for small 
propagation distances, is explained here by the fact that off-diagonal elements of H are 
unmodified by matter. The result, Eq. ^6], is valid for all mass and mixing patterns, 
but is, of course, useful only for small propagation distances. 

A particular success of the matter invariance, Eq. [IB], is that it enables the calcu- 
lation of the universal CP violating asymmetry in matter of uniform density, with no 
restriction to small L, knowing only the vacuum Hamiltonian and the matter-modified 
eigenvalue differences (which are readily calculated in terms of vacuum quantities, 



1 ^ 1 ^ ' (A' 12 /2)(A 23 /2)(A' 31 /2) 

(27) 

This expression does not require knowledge of the matter-modified eigenstates (ie. mix- 
ing parameters), whose dependence is contained completely within the matter- invariant 
factor Im(P eM P AtT .P Te ). A comparison of the CP-conserving and CP- violating terms 
in Eq. ^ leads us to suspect that, perhaps more generally, the matter-invariant quan- 
tities \H a/ 3\ 2 play a similar role for CP-conserving observables to that played by 
\m{H eil H^ T H Te ) for CP-violating observables. It is interesting therefore to enquire 
whether a simplification similar to that of Eq. |27| exists also for the CP-conserving 
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parts of P(a — > (3), in terms of the matter-invariants \H a p\ 

Some guidance is obtained by returning briefly to the two-family case. Consider 
Eq. |] for the e — > \i appearance probability in matter of uniform density: 

P'(e-/i) = |^/ = -4<Jfsin 2 (A 21 L/2) 

sin 2 (A' 21 L/2) 
" (A 21 /2) 2 (28) 

where we have used Eq. ^ to write the probability as the product of a matter-invariant 
part (\H' efM \ 2 = \H e ^\ 2 ) and a part with a generic dependence on eigenvalue differences 
of the form [sin (A'L/2)/(A'/2)] n . This exact expression, Eq. |28|, is clearly analogous 
to Eq. ^7] in the sense we require. 

We do not a priori expect the three-family case to be as straightforward, but it 
will turn out that much of the above simplicity is retained, especially in the case of 
a hierarchical neutrino spectrum (as seems to be preferred experimentally). We may 
start again with Eq. |8], and assume a hierarchical spectrum, A' 21 << A 31 , A 31 ~ A 32 = 
A' for all matter-densities considered. Then in the limit A' 21 L/2 << 1 (which is less 
restrictive than the limit in which Eq. |26] applies), we find for the leading oscillation of 
the appearance probability (disappearance probabilities are anyway readily obtained 
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from appearance probabilities through unitarity): 



- 0) = ~ -4«5 + < 3 |) sin 2 (A'L/2) 

= 4|x2 5 /3 | 2 sin 2 (A , L/2) (29) 

where the last step follows from unitarity. We also require that jXAlA^ << |X^|A' 31 
for all matter-densities considered, in order that the leading oscillation dominates. This 
will generally be true, for the hierarchical spectrum we have assumed, but may not be 
in some special cases eg. if the mixing matrix elements also have an extreme hierarchy, 
see below. 

Applying the same approximations to Eq. |21| we have: 
\X'i\ 2 A- = + <g)A' 2 * \H' ap \* = \H aP \* * -(Kll + i^)A 2 = |Xi| 2 A 2 (30) 

(ie. the combination |X^|A' is a matter-invariant, as long as (X^A^ << (X^A^ 
remains valid, as already assumed). Combining Eqs. ^ and |3D], we find that 

V>(„ R\ \A> I* -Iff ,2 sm 2 (A'L/2) 

P{a-+(3) = \A a/3 \ ~ \H afS \ , A ,i 2 \2 ( 31 ) 

which generalises Eq. ^8| to the three-family case, when the specified approximations 
apply. 



Thus Eq. |31] comes close (as close as we believe is possible) to being the analogue 
of Eq. |27] for the CP-conserving case. Where valid, Eq. |3l| retains all the calculational 
expediency previously recognised in connection with Eq. |27|, for the CP-violating case. 
We should say that in general (for arbitrary mass-squared differences, mixing angles 
and propagation lengths) no such analogy exists between the CP-conserving and the 



CP-violating contributions. The simplification exemplified by Eq. |2j is a consequence 
of the universality (for three families only) of the CP-violating term, as embodied 
in the Jarlskog parameter, J ||, but in general, no such universality exists for the 
CP-conserving K 1 ^ coefficients (the two-family case being an exception). However, 
for phenomenologically realistic oscillation parameters (see eg. [0), Eq. [IT] is cer- 
tainly applicable to the leading oscillation of P'(/i — > r) (and of P'(r — > //)) for most 
energies of interest, and, unless it turns out that |JfgJA 3 i ~ jX^JAai for (3 = fi, r 
(|-X"efll depends on the presently poorly determined U e 3), it should also be relevant for 
appearance probabilities involving electron neutrinos. 



4 Relation to Standard Formalism 

In order to relate our matter-covariant formulation to the standard one in terms of 
mixing angles, we first need to remove the implicit dependence of our matter-invariants 



3 The (possibly realised |13|]) case |f/ e 3| 2 = 0, is a specific case which violates the assumption 
l-XjjglAju << jXjglA^ for a = e or (3 — e, because then, X\ — 0, and the leading oscillation 



vanishes. However, in this case, the long-wavelength oscillation becomes dominant, and Eqs. 29 and 

1 — ¥ \X a p\ 



become exactly true, with the substitutions \XjL\ — > \X^g\ and A' = A^ etc. 
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on the neutrino energy (which has been left in until now in the interests of a concise 
discussion). We note that all the preceding equations in this paper remain unaltered, 
if we make the energy re-scaling H — > 2EH = TC, A^ — > 2EAij = Am?-, \ — > 2EAt, 
a — > 2Ea = A, L —>■ L/2E and likewise for all primed quantities. The observables 
K^p, J, and their matter- modified counterparts remain unchanged throughout, as do 
the oscillation probabilities. 

We can now give the relationships between our energy-scaled matter-invariants 
(and covariant) (Eqs. and [T9| - f23"|) and the six parameters conventionally used to 
describe neutrino oscillation phenomenology, #12, 623, #13, 5, Amjj and Am^. Defining 
D = Am^ — Am^Sy and F = ^Si 2 S23Si 3 Ci 3 , we have: 

Hi = \H e ,\ 2 = c 2 13 [s 2 23 s 2 13 D 2 + \c 2 23 S 2 12 {Am 2 21 ) 2 } + 2FDAm 2 21 c 5 (32) 

H 2 = \H eT \ 2 = c 2 l3 [c 2 23 s\ 3 D 2 + \s 2 23 S 2 12 {Am 2 21 ) 2 } - 2FDAm 2 21 c s (33) 

H 3 = \H, T \ 2 = ~{[c 2 l3 S 23 D - (C 12 S 2 3 + S l2 C2 3 s 13 cs)Am 2 21 ] 2 + S 2 12 s 2 13 s 2 (Am 2 21 ) 2 } (34) 

Hi = lm(n e ^n^ T n T e) = Am 2 21 Am 2 31 (Am 2 31 - Am 2 21 )Fs 5 (35) 
^5 = - H TT = C 2 2,[C 12 Aml 1 - c\ 3 D] - S 12 S 2 3Si3C S Aml 1 (36) 

Ho = H ee - l -TrH = D(s 2 l3 - ~) - l -C 12 Am 2 2l (37) 



where we have used the conventional parameterisation |14| and c^- = cos^-, Sy = sin^-, 
Cij = cos28ij, Sij = sm28ij etc. Eqs. P^-P^| remain true in matter (with primed 
quantities on the RHS), as well as in vacuum and show that (sums of) products of 
(matter-density-dependent) oscillation amplitudes and mass-squared-differences are 



matter- invariant, thereby accounting for the folklore "no-win" theorem ||15|| , that 
increased mixing-angles in matter seem to be compensated by reduced masses and 
hence longer oscillation lengths, and vice-versa. The covariant quantity on the LHS 
of Eq. |37| transforms to the matter H'q = K' ee - Tr(H')/3 = H + 2A/3, 

where A = 2\^2GpN e E, while the quantities on the RHS take their density- dependent 
values. 

The usual six parameters (6*12, 9 23 , #13, 5, Am^ and Am^-J, are popular for their 
intuitive relevance, but they all have complicated dependencies on the density of the 
matter through which the neutrinos propagate [0 [|T^] . Different experiments operating 
over different baselines at different energies are all subject to different modifications 
of each of these conventional parameters due to matter effects. However, in each case, 
only one additional parameter is actually present, the (path- averaged) Wolfenstein 



term, A. By contrast, as we have seen, the five combinations of observables, Eqs. |32- 
|36| are all matter-invariant, while the sixth observable, Eq. |37], has a very simple 
dependence on the matter density. Clearly, these quantities could be used as the basis 
for a new parameterisation of neutrino oscillation phenomenology. This is especially 
true for the matter- invariants involving only off-diagonal elements: \Ti. a f3\ 2 (a 7^ (3) 
and Im(7i e/J 7i /JT 7i re ), which are determined directly in appearance experiments as the 
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amplitudes of the leading CP-even oscillations (given a strong-enough mass-hierarchy) 
and the CP-asymmetry respectively. Not only are they matter- invariant, but they are 
more directly observed even than the conventional mixing angles. 

To some degree then, one may say that the problem of deconvoluting matter effects 
on the standard mixing parameter set is a problem of man's own making, which is 
solved (at least in the uniform-density case we have considered) by simply switching 
to a different parameter set. While it is probably unrealistic to assume that the new 
parameterisation proposed here will be immediately adopted universally, we believe 
that our parameters may have a practical value, beyond their evident conceptual 
importance as matter-invariant observables. 
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